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Abstract 

The cross section of the e + e~ — > ujtt — > 7r°7r°7 reaction was mea- 
sured by the SND detector at VEPP-2M e + e~ collider in the energy 
range from threshold up to 1.4 GeV. Results of the cross section fitting 
by the sum of p, p' and p" contributions are presented. 



1 Introduction 

The process of e + e~ annihilation into hadrons in the 1-2 GeV energy region 
is an important source of information about excited states of light vector 
mesons p, uo and 0. Parameters of these states are still not well established 
mainly due to poor accuracy of existing experimental data. Recently new 
data in this energy region became available. The cross sections of the re- 
actions e + e _ — > 37r Jl| and e + e~ — > An were measured at VEPP-2M 
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collider. Accuracy was also significantly improved in the recent measure- 
ments of r —>■ 2tti> t and r — ► Attu t spectral functions || |4], |5| related to 
corresponding e + e~ annihilation cross sections by the CVC hypothesis W\ . 

The reaction e + e~ — > urn is one of the dominant processes in the energy 
range from 1 to 2 GeV. The PDG values for the mass and width of the 
p'(1450) meson are based mainly on phenomenological studies of this 
process ||. The most precise measurements of the e + e~ — > ujtt cross section 
and t — > u;7rz/ r spectral function were done in 0, §, f4j. All these experiments 
addressed the 4ir final state into which other intermediate states, for example 
ai7r, contribute significantly. In this work the e + e~ — > co>7r° reaction was 
studied in the 7r°7r°7 final state where other contributions are much smaller. 
This allows to avoid systematic errors inherent to the e + e~ — > 47r channel 
due to non-trivial background subtraction, which must take into account 
interference effects. The process e + e~ — ► co>7r° — ► 7r°7r°7 was studied for the 
first time in M where 20% statistical accuracy was achieved. 



2 Detector and experiment 

SND is a general purpose non-magnetic detector [|Kj . Its main part is a three- 
layer scintillation electromagnetic calorimeter consisted of 1630 Nal(Tl) crys- 
tals with solid angle coverage about 90% of 47r. The energy resolution of the 



calorimeter for photons is <te/E = 4.2%/ <g E(GeV), the angular resolution is 
about 1.5°. The directions of charged particles are measured by two coaxial 
cylindrical drift chambers covering 95% of 47r solid angle. 

The analysis presented in this work is based on data recorded in 1997- 
1999 in two separate energy regions: 920-980 MeV and 1040-1380 MeV. 
Analysis of the ^-resonance region (980-1040 MeV) was published earlier 
p4| . In the first energy region the total integrated luminosity of 1.5 pb _1 
was collected at 6 energy points. The region above the <fi meson was scanned 
with a 10 MeV step. Total integrated luminosity accumulated in this region 
is about 9 pb _ . The luminosity was measured with a systematic uncertainty 
of 3% using e + e~ — > e + e~ and e + e _ — ► 77 reactions. 

3 Event selection 

For primary selection of 

e+e- -> CJ7T -> 7r°7r°7 (1) 
events the following criteria were applied: 
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Figure 1: The Xniry distribution for 
the experimental (points with error 
bars) and simulated (the histogram) 
events of the process ([[]). 



Figure 2: The distribution for 

the experimental (points with error 
bars) and simulated (the histogram) 
events of the process ([[]). 



• five or more photons and no charged tracks are found in an event; 

• the energy deposition in the calorimeter is more than 0.7.E; 

• the total momentum of an event measured by the calorimeter is less 
than 0.15E; 

where E is a center of mass energy of e + e~ pair. The main sources of back- 
ground surviving these cuts are QED processes 

e + e- -> 2 7 , 37. (2) 

with extra photons either from the the beam background or splitting of elec- 
tromagnetic showers in the calorimeter. 

For each event satisfying primary selection criteria the kinematic fitting 
assuming e + e _ — > 7r°7r°7 — > 57 hypothesis was performed. As a result, two 
parameters: Xmry — the x 2 °f the fit and M ny — the 7r° recoil mass closest 
to that of to meson were evaluated. 

The x-K-K-y distribution for the events with \M^ — 782 1 < 50 and the n° 
recoil mass spectrum for the events with x-mry < 30 are plotted in figs.[j] 
and |2|. Good agreement between experimental distributions and simulation 
of the process (0) shows that there are no other significant contributions in 
the selected event sample. For final event selection the following cuts were 
applied: 

X™ 7 < 30 > l M -7 - 782 l < 50 > ( 3 ) 
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which reject the QED background (0) almost completely. In order to estimate 
residual background the experimental 7r° recoil mass spectrum was fitted by 
the sum of simulated spectrum of the process ([[]) and a linear background. 
As a result the estimated total number of background events did not exceed 
1.5% of all selected events. This value was taken as a systematic error of the 
measured cross section of the process (|l|) related to the residual background. 

Experimental data were collected in 40 energy points. In the energy region 
under study the cross section changes slowly so it was possible to reduce 
the number of energy points combining the neighboring ones. Resulting 
energies, their standard deviations, integrated luminosities (IL) and numbers 
of selected events (N) are listed in the Table [I]. The detection efficiency for 
the process (|T|) was determined by simulation based on a formula from \T2 
which takes into account finite width of a>-meson. In the energy range from 
1050 to 1400 MeV the efficiency was found to be constant and equal to 
40%. Near the e + e~ — > um° threshold the number of events with the tt° 
recoil mass below cu-meson mass increases sharply and the efficiency goes 
down. The energy dependence of the detection efficiency is presented in the 
Table 0. 



4 Fitting of the cross section 

The visible cross section a V i S = N/IL is related with the Born cross section 
of e + e~ — > um° — > 7T 7r°7 process as 

a vis (E)=e(E)a (E)(l + 5(E)), (4) 



where 5(E) is a radiative correction calculated according to [13] . Radiative 
corrections for the different energies and obtained cross section values are 
listed in the Table [I]. The cross section at E — was taken from [11 



Only statistical errors are shown in the table. The systematic error includes 
the error of the luminosity measurement (3%), the detection efficiency error 
(4%), possible background contribution (1.5%), and the error of radiative 
correction (1%). The total systematic error was estimated to be 5%. 

Our results in comparison with the most precise CMD-2 ||, CLEO 
and DM2 || measurements are shown in Fig|| The cross sections from 
II U measured in the 7T + 7r~7r°7r channel were recalculated using the PDG 
value of B(u) — > tt 7)[0]. The cross section from |4j] was obtained from the 



r — > uj-KVr spectral function assuming the CVC hypothesis |14| . The CLEO 
results are in good agreement with ours while the CMD-2 measurements are 
about 10% lower, although the difference observed is smaller than the 15% 
systematic error quoted in H. There is a significant difference between the 
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Table 1: Energy, its standard deviation, integrated luminosity, detection efficiency, 
number of events, radiative correction and cross section of the process e + e _ — ► 



E, MeV 


AE, MeV 


IL, nb 


e 


N 


8 


(To, nb 


920 


0.2 


328 


0.169 


3 


-0.148 


0.06±0.04 


940 


0.2 


289 


0.312 


15 


-0.168 


0.20±0.05 


954 


2.0 


496 


0.332 


39 


-0.157 


0.28±0.05 


973 


4.3 


373 


0.352 


68 


-0.140 


0.60±0.07 


1020 


— 


- 


— 


— 


- 


0.74±0.02 


1045 


5.0 


152 


0.390 


51 


-0.100 


0.96±0.13 


1063 


4.4 


371 


0.390 


124 


-0.093 


0.95±0.09 


1081 


3.5 


666 


0.390 


253 


-0.086 


1.07±0.07 


1102 


3.8 


524 


0.390 


180 


-0.080 


0.96±0.07 


1123 


4.4 


420 


0.390 


171 


-0.074 


1.13±0.09 


1142 


4.0 


358 


0.390 


142 


-0.070 


1.09±0.09 


1160 


0.4 


316 


0.390 


158 


-0.066 


1.37±0.11 


1183 


4.6 


587 


0.390 


280 


-0.061 


1.30±0.08 


1203 


4.4 


569 


0.390 


284 


-0.057 


1.36±0.08 


1223 


4.6 


465 


0.390 


232 


-0.054 


1.35±0.09 


1243 


4.8 


562 


0.390 


309 


-0.050 


1.48±0.08 


1266 


5.1 


397 


0.390 


241 


-0.046 


1.63±0.11 


1285 


5.2 


492 


0.390 


268 


-0.043 


1.46±0.09 


1304 


5.0 


459 


0.390 


249 


-0.040 


1.45±0.09 


1325 


5.0 


516 


0.390 


301 


-0.038 


1.55±0.04 


1344 


5.8 


676 


0.390 


387 


-0.036 


1.52±0.08 


1363 


4.6 


857 


0.390 


497 


-0.034 


1.54±0.07 


1380 


0.5 


470 


0.390 


234 


-0.033 


1.32±0.09 



results of DM2 |J and CLEO ||. For the cross section fitting we used our 
data together with the data from CLEO. 

The energy dependence of the process (H) cross section was written as a 
sum of contributions from p(770) and its excitations p' and p": 



a {E) 



Ana 2 (gpun 



E 3 V /, 



rn £Lr 



m„,. 



D r , 



p"U)-K 



P f {E). (5) 



Here a is a fine structure constant and g puJ7T is a p — > ujtt coupling con- 
stant. The f p coupling constant was calculated from the p — > e + e~ decay 
width: r pee = 4nm p a 2 /3fp. The expression m 2 /D p represents p-meson Breit- 



Wigner amplitude with D, 



mz 



■E 2 —iET p (E), where m p and T p (E) denote 



the p meson mass and energy-dependent total width respectively. The real 
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parameters A\ = g Pi(M / 'g^-n • fp/fpi are the ratios of the coupling constants 
of different p states. The factor Pf(E) describes the energy dependence of 
the final state phase space. In the case of infinitely narrow u resonance 
Pf(E) = 1/3 ■ q^B(uj — > 7r°7), where q u is an cj-meson momentum. This 
approximation is good for the energy range above 1050 MeV, but at energies 
close to the e + e~ — > ujtc threshold it is more adequate to use precise formula 
taking into account finite width of the u meson ||12|| . The Blatt-Weisskopf 
factors C PiUJn restricting fast growth of the partial widths, were taken in the 
form IBfl: 



1 + {Rq.{E)r ^ 1 + (Rq.(E)) 2 



The range parameter R was assumed to be the same for p, p' and p" mesons. 
The energy dependence of the p(770) total width was expressed as: 



UE) - r pK) ^j {^)C%« + ^E)CU (7) 

where q n is the pion momentum, C pn7T is the Blatt-Weisskopf factor for p 
7t + tt~ decay |15| : 



C p7T7T (E) 



1 + {Rq,{m p )Y 



\ 1 + (RqAE)) 2 ' 

There is no generally accepted description of the shapes of the broad excited 



states p' and p" . In [|] constant total width were assumed, while in Jl6[ the 
total widths varied with energy as a sum of the partial widths of all main 
decay modes. We considered the both approaches in order to understand 
model dependence of the fit parameters. 

The fit results obtained for three classes of models are listed in the Table 0. 
The fit parameters are sensitive to the variation of the range parameters R 
in the Blatt-Weisskopf factors. In the Table we show the intervals of the fit 
parameter variation when R ranges from to 2 GeV -1 . The typical errors 
of the parameters obtained for each specific model are listed in the fourth 
column of the Table 

In the model 1 the constant widths of p' and p" were assumed. The 
mass and width of the p" meson were fixed to their PDG values [0: m p n = 
1700 MeV, Tp" = 235 MeV. The obtained value of the p" amplitude A 2 
is compatible with zero, so the final fit result for the model 1 is given for 
A 2 fixed to zero. Let us discuss the choice of an upper boundary for the R 
parameter. Small values g puJ1T < 13 GeV -1 found for R > 2 GeV -1 are in 
conflict with the QCD sum rules estimation 16 GeV -1 [17] and experimental 
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Table 2: The fit results for various models described in the text. 





Model 1 


Model 2 


Model 3 


Error 


R, GeV -1 


0-2 


0-2 


0-2 




9pum, GeV" 1 


16.1-13.3 


15.6-13.2 


15.6-12.9 


0.3-0.4 


Try, MeV 


1460-1520 




= 1400. 


10 


T//, MeV 


380-500 




= 500. 


20-30 


Ai 


-(0.22-0.24) 


= 0. 


-(0.04-0.08) 


0.01-0.02 


m p //, MeV 




1710-1580 


1620-1550 


15-20 


I>, MeV 




1040-490 


580-350 


70-20 


A 2 


= 0. 


-(0.22-0.20) 


-(0.18-0.13) 


0.1 


x 2 /nd 


(52-48)/35 


(47-48)/35 


(43-44)/34 





value of 14.4 GeV -1 , obtained from the u — > 3tt decay width assuming that 
uj — > pir mechanism dominates in this decay. Only for R ~ the extracted 
mass and width of the p' meson are compatible with their PDG values 0. 
However in this case the fit yields the lowest confidence level of P(x 2 ) — 3%. 



Models 2 and 3 take into account energy dependence of the total widths 
of p' and p" mesons. Since branching ratios of p' and p" decays are practi- 
cally unknown @ the energy dependence of the sum of all rapidly growing 
partial widths of their multihadron decays was approximated as the energy 
dependence of pi — > ujtx width: 



l-B 



+ B. 



ft— >7T7T 



m. 



E 



C 



PiUITT 



r 2 

Pi-KTT 



(9) 



For the p" 
10% 



— > 27r branching fraction we use the theoretical estimation: B p /r^ n7T ■■ 
For the p' meson the value B p /^ n7T = 50% was chosen which, we 
think, reflects experimental situation more correctly than the theoretical pre- 
diction 25% |T3 . The only difference between the models 1 and 2 which both 
consider only one excited p state is that the model 2 assumes its total width 
energy dependent. But the fit results for these two models are quite differ- 
ent. Particularly the mass found for the model 2 is close to the PDG value 
of the p"-meson mass. On the other hand there is a definite signal of the p' 
meson with a mass of 1320-1400 MeV and 400-500 MeV width in the pion 
form factor data || 0. Therefore we also considered the model 3 with two 
excited p states, in which the mass and width of the p' meson were fixed 
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1000 1200 1400 1600 1800 2000 2200 

E, MeV 

Figure 3: The cross section of the reaction e + e~ — > wzr — > 7r°7r°7. The results of 
the SND (this work), DM2 §§, CMD @, CLEO § experiments are shown. Curves 
are results of fitting to the data in model 1 and 3 with R = 0. 



to m p i = 1400 MeV and I> = 500 MeV. This model yields the best fit to 
the experimental data: P{x 2 ) ~ 13%. Fitting curves corresponding to the 
models 1 and 3 with R = are shown in figure [| 

The main conclusions from the analysis of the fit results are the following. 
Fitting of the same experimental data by models with fixed and energy- 
dependent total widths of the excited states yields quite different parameters 
of these states. This is caused by strong energy dependence of the phase space 
for the main decay modes of p' and p" mesons and this effect should be taken 
into account in the fitting of experimental data. Satisfactory description of 
the experimental cross section was obtained in the model with two excited 
states with the masses m p i = 1400 MeV and m p » ~ 1600 MeV in which 
contribution of the higher state dominates. However this result contradicts 



the theoretical expectation |1| pLgf] , where p' and p" are considered as 25" 
and ID qq states respectively and the larger contribution of the lower 2S 
excitation was predicted. Thus, with the new experimental data the situation 
in the isovector sector remains unclear. The main problem for data analysis 
is the absence of consistent phenomenological description of the shapes of 
broad resonances with strong energy dependence of partial widths. 
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5 Summary 



In this work the cross section of the e + e~ — > ujtt — > 7r°7r°7 reaction was 
measured from the threshold up to 1.4 GeV with a 5% systematic accuracy. 
This is the most precise measurement of the e + e~ — > kJ7r° cross section in this 
energy range. Our data are in a good agreement with the CLEO measurement 
of t — > co>7rz/ T spectral function |3J]. The combined fit to our and CLEO data 
was performed in the vector meson dominance model taking into account the 
contributions of p, p' and p" states. The experimental data can be reasonably 
well approximated assuming existence of only one excited state but its mass 
is strongly model-dependent and varies from 1460 to 1700 MeV for different 
descriptions of the resonance shape. The best fit to the data was obtained for 
the model with two excited states in which the higher state with m p » pa 1600 
MeV dominates. 

This work is supported in part by the Russian Fund for Basic Researches 
(grants No. 99-02-16815, 99-02-17155) and STP "Integration" (grant No. 
274). 
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